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A technique for detecting brain injury at the bedside has great cUnical value, but conventional imaging 
techniques (such as computed tomography [CT] and magnetic resonance imaging) are impractical. In this 
study, a novel method-the symmetrical channel electroencephalogram (EEG) signal analysis-was developed 
for this purpose. The study population consisted of 45 traumatic brain injury patients and 10 healthy 
controls. EEG signals in resting and stimulus states were acquired, and approximate entropy (ApEn) and 
slow- wave coefficient were extracted to calculate the ratio values of ApEn and SWC for injured and 
uninjured areas. Statistical analyses showed that the ratio values for both ApEn and SWC between injured 
and uninjured brain areas differed significantly (P < 0.05) for both resting and name call stimulus states. A 
set of criteria (range of ratio values) to determine whether a brain area is injured or uninjured was proposed 
and its reliability was verified by statistical analyses and CT images. 

Quantitative and topographic analysis of the human electroencephalogram (EEG) is a valuable method for 
evaluating brain function in healthy and diseased states, as those analyses could permit determination of 
scalp signal sources. For example, it has proven useful methods for mapping the seizure patient's 
"epileptogenic zone"^'^. It could also be used to characterize different brain pathological states, such as psychiatric 
and neurological disorders^. However, whether quantitative and spatial EEG analysis can be used as a technique 
for mapping injured brain areas has not been sufficiently explored^; recent studies that have attempted to resolve 
this issue^'^, and the current investigation is based on the conclusions from these reports. 

Injured regions in the brain are typically identified using computed tomography (CT) or magnetic resonance 
imaging (MRI). However, compared with CT or MRI, the EEG method has advantages for its cost, high temporal 
resolution, flexibility, ease-of-use, performance, and speed As such, using EEG at the bedside can be a powerful 
tool that provides dynamic information about brain functions^"^, and if applied to the assessment of brain injury, 
has great theoretical and clinical value^°'^\ 

In this study, an EEG analysis method was developed that is refered to as symmetrical channel EEG signal 
analysis (SESA). The SESA method is based on hypothetical differences in EEG signal features between injured 
and uninjured brain areas (applicable for unilateral brain injuries, the most common brain injuries). If the 
difference is statistically significant, it could be used as an index to judge whether or not a brain area is injured. 
The SESA method is based on the characteristic anatomical and neurophysiological symmetry and functional 
offside symmetry of the brain 

Previous studies have shown that EEG signals for corresponding regions in the cerebral hemispheres have a 
high degree of similarity in the resting state^^'^^. If a unilateral brain injury occurs, we can infer that this EEG 
similarity will decrease and differences will increase in the symmetrical cerebral hemisphere. In the active state, 
characteristic offside symmetry also increases EEG signal differences^^. Indeed, recent developments in brain 
computer interfaces suggest that pattern recognition of imagined movements are based on the characteristic 
functional offside symmetry exhibited by the brain^^"^^. Therefore, it is feasible, both in theory and practice, to 
map brain injury areas using the SESA method. 

We applied signal processing methods and statistical analyses to verify whether changes in EEG signal features 
(in the symmetrical cerebral hemisphere) are real, and whether they can be used as a method for determining 
brain injury. A set of criteria (range of ratio values) to distinguish between injured and uninjured brain areas was 
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Table 1 1 
Features 


Mean Cp values for total injured 
oTares 


and uninjured channels 

Mean Cp tor injured chiannels 


Mean Cp tor uninjured chiannels 


P value 


ApEn 


Resting 


0.665 ± 0.007 


0.905 ± 0.009 


1.98E-13** 


Name call 


0.584 ± 0.006 


0.874 ± 0.005 


3.05E-12** 


SWC 


Resting 


1.711 ±0.047 


1.343 ± 0.017 


5.25E-07**s 




Name call 


1.957 ± 0.166 


1.517 ±0.072 


9.8E-06** 


One-way analysis of variance (ANOVA) for control groups showed Rvalues all > 0.05 (not listed in the table), indicating that there were no differences between symmetrica 
for control and n = 45 for patient groups, respectively. 

**P < 0.01 for injured vs. uninjured channels. ApEn, approximate entropy; Cp, ratio of signals from injured and uninjured areas; SWC, slow-wave coefficient. 


channels for the controls, n = 1 0 



also proposed based on SESA. To test whether the proposed criteria 
were rehable, statistical vaHdation was carried out, and the corres- 
pondence between diagnoses made based on SESA criteria and CT 
images was compared for five selected patients. 

Results 

The SESA method was developed as follows. First, two EEG signal 
features-approximate entropy (ApEn) and slow-wave coefficient 
(SWC) were extracted for all brain channels, and the ratio (Cp value) 
of signals from injured and uninjured areas for these features was 
calculated. Statistical analyses of Cp values in population and indi- 
vidual channels were carried out, and a range of ratio values was set as 
a standard for determining whether a brain area is injured or 
uninjured. Finally, the reliability of the SESA method was validated 
using CT images. 

Statistical analysis of Cp values for total injured and uninjured 
channels. A statistical analysis was performed to evaluate differences 
in population Cp values for total injured and uninjured channels. 
The results are shown in Table 1. 
Results from Table 1 indicated that: 

(1) For ApEn and SWC, in both the resting and active states, a P < 
0.01 was observed for population mean Cp values, indicating 
that injured and uninjured brain channels differed signifi- 
cantly; 



(2) For the Cp value of ApEn, the value of the injured area during 
the call name stimulus was much smaller than the value in the 
resting state; 

(3) For the Cp value of SWC, the value of the injured area during 
the call name stimulus was much larger than the value in the 
resting state. 

The Cp value for ApEn during the name call stimulus was smaller 
than in the resting state. As described in the Methods, in the resting 
state, the ApEn values for both the injured (e.g., Fpl) and symmet- 
rical uninjured (e.g., Fp2) channels did not change significantly (i.e., 
Fpl/Fp2 remained constant); however, during the name call stimu- 
lus, the ApEn value for Fp2 increased but was unaltered for Fpl due 
to injury, resulting in an overall decrease in Fpl/Fp2. Conversely, the 
Cp value for SWC during the name call stimulus was larger than in 
the resting state because Fp2 was unchanged while Fpl increased. 

Statistical analysis of Cp values for individual injured channels. 

An analysis of population Cp values for all injured and uninjured 
channels revealed differences in the Cp values for individual injured 
channels (Figure 1). 

The results showed that: 

( 1 ) For ApEn and SWC, in both resting and active states, there were 
statistically significant differences (P < 0.01) for all individual 
injured channels relative to uninjured channels; 
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Figure 1 | Cp values for individual injured channels. Values for ApEn are shown; those for SWC can be found in Appendix !.**?< 0.05. 
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Table 2 


Range of Cp values for injured channels 










Injured Channel 




Features 


otare rp i or rpz 


r/or rO 1 o or 14 


1 0 or 10 


ApEn 
SWC 


Resting 0.583 ~ 0.699 
Call Name 0.503 ~ 0.638 
Resting 1 .554 ~ 1 .952 
Call Name 1.736-2.466 


0.626 ~ 0.744 0.593 ~ 0.753 
0.520-0.638 0.543 -0.665 
1.615 ~ 1.935 1.537- 1.778 
1.612 -2.161 1.587-2.131 


0.629-0.717 
0.535 -0.636 
1 .502 - 1 .694 
1.755 -2.164 


A brain area was considered as injured if the corresponding Cp values were within the ranges indicated for F3, F4, P3, P4, C3, C4, Ol , and 02 channels. ApEn, 
coefficient. 


approximate entropy; SWC, slow-wave 




Table 3 


Range of Cp values for uninjured channels 










Uninjured Channel 




Features 


State F3 or F4 


C3orC4 P3orP4 


Ol or02 


ApEn 
SWC 


Resting 0.634 ~ 1 .071 
Call Name 0.627- 1.125 
Resting 0.997 ~ 1 .394 
Call Name 1 .056 ~ 1 .545 


0.722 - 1.106 0.737- 1.165 
0.697- 1.127 0.736 - 1.216 
1.001 - 1.445 1.026- 1.453 
1.1 1 1 - 1.709 1.146 - 1.61 1 


0.716- 1.233 
0.650- 1.413 
0.952 - 1 .474 
1.1 18 - 1.651 


A brain area was considered as uninjured if the corresponding Cp values were within the ranges 
coefficient. 


indicated for F3, F4, P3, P4, C3, C4, Ol , and 02 channels. ApEn 


approximate entropy; SWC, slow-wave 



(2) The Cp value of ApEn for injured channels during the name call 
stimulus was smaller than in the resting state; 

(3) The Cp value of SWC for injured channels during the name call 
stimulus was larger than in the resting state (details in 
Appendix III); 

An explanation for why the Cp value decreases for ApEn while 
increasing for SWC is given above (see Appendix I). 

Criteria for assessing injured vs. uninjured brain channels. The 

above results indicated that there are significant differences in Cp 
values between injured and uninjured areas as well as individual 
channels. Thus, statistically determined differences served as a 
basis for establishing criteria that distinguished between injured 
and uninjured brain areas. These criteria (section 7 of Methods) 
and preliminary ranges of Cp values are given in Tables 2 and 3. 

Verification of accuracy. Permutation testing was carried out to 
assess the accuracy of the established method. For example, the Cp 
values of ApEn for Fpl or Fp2 with the name call stimulus comprised 

0.96 



31 cases; of these, 20 were randomly selected as a training group, and 
the remaining 1 1 were the test group. The detailed procedure was as 
follows: the range of Cp values for the training and test groups were 
calculated using the above- described method; the reliability was 
assessed by determining the region of overlap between the Cp 
values of the training and test groups. This process was repeated 
for randomly selected training and test cases. After 50 iterations, 
the average accuracy was found to be 92.3%. 

The calculation was repeated for the injured leads of F7 or F8, T3 
or T4, and T5 or T6, and uninjured leads P3 or P4, C3 or C3, F3 or F4, 
and Ol or 02 (Figure 2). The results showed a high average accuracy, 
confirming that the method is reliable. 

Calculation of probabilities of type I and II errors. The probability 
of occurrence of type I and II errors (incorrectly judging an uninjured 
area as injured one, or vice versa) was determined. The Cp values of 
ApEn for Fpl or Fp2 is used as an example. Of 45 patients, 31 were 
injured at Fpl or Fp2, and 14 were uninjured in these locations. 
Assuming that the 14 cases follow a normal distribution, Cp values 
of ApEn for Fpl or Fp2 were calculated as for injured cases. The 



ApKn 
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Figure 2 | Verification of criteria for distinguishing between injured and uninjured brain areas. Data represent the average accuracy after 50 iterations of 
permutation testing. The lowest average accuracy was from the T5/T6 leads, partly due to a smaller sample size for this group (12 cases). 
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probability 




Probability 
of type II 
error 



Probability 
of type I 
error 



Cp value 



Figure 3 | Method for calculating the probabilities of type I and II errors. 

The probability distribution of ApEn Cp values for injured and uninjured 
cases is shown for the Fp 1 or Fp2 leads. 

distribution of the ApEn Cp values for the two groups is shown in 
Figure 3. The probability of type I and II errors was calculated at the 
point at which the curves intersect and determined as the areas under 
the curve. 

The calculated probabilities for all leads are shown in Table 4. 
Since the values were small, it can be concluded that the method is 
reliable. 

Verification of the SESA method using CT images. To test the 
reliability of the proposed criteria, five additional patients were 
selected for whom diagnoses made based on the SESA method and 
CT images were compared. The Cp values from one patient are listed 
in Table 5. In this case, the exact location of the injured area is 
unknown. Calculating Cp values of the symmetric channel 
permitted determination of the injured areas. According to the 
diagnostic criteria, brain areas corresponding to channels with Cp 
values shown in bold italics in Table 5 (FP1/FP2 and F7/F8) were 
injured. 

An example of results obtained using the SESA method and CT 
images for one patient is shown in Figure 4. A good correspondence 
was observed between the two methods. 

Discussion 

The SESA method was developed to map injured brain areas, and its 
clinical utility was verified using theoretical and experimental meth- 
ods. Although theoretical simulation studies have been conducted in 
which functional analyses of asymmetric EEC recordings have been 
applied towards clinical assessment^^"^^, our report investigates in 
practice the feasibility of using EEC signals for functional brain 
mapping. As an innovative brain mapping method, it has the poten- 
tial to be used at the bedside in the clinic. Through further develop- 
ment, it could also serve as a powerful approach or portable imaging 
device for brain function analysis. Other fundamental work, includ- 
ing more analysis tools and more engineering, must be carried out to 
develop its practical use. We will specifically discuss some fun- 
damental work required. 

The scope of SESA applications. Brain injuries can have different 
causes, such as cerebral trauma and cerebral ischemia. From a 



topographical perspective, brain injuries can be classified as 
unilateral injuries, bilateral injuries, dispersed injuries, and whole 
brain injuries. The SESA method is especially suitable for 
unilateral brain injuries, as it is based on hypothetical EEC signal 
differences between asymmetrical signals of brain injury and 
symmetrical signals for uninjured areas. 

In the clinic most brain injuries are unilateral; however, whether 
SESA is suitable for analyzing other kinds of brain injuries requires 
further investigation. To this end, a method for multi-area neural 
mass modeling based on EEC and MEG signals has been developed^^. 
We suggest that other methods, such as regional symmetrical EEC 
signal analysis, might also have the potential for applications to 
unilateral and whole brain injuries. 

In addition, our research sample included patients with severe 
disorders of consciousness. Compared with brain injury patients 
with normal consciousness, these patients' EEC signals are relatively 
stable^^'^^, and thus used for this study. Further research will be 
performed to experimentally test and verify whether this method is 
suitable for assessing brain injuries associated with normal 
consciousness. 

Establishment of diagnostic criteria using patient data. A range of 
ratio values defined for each channel is proposed as a standard for 
determining whether or not a brain area is injured. Although the 
method was validated by statistical analyses and was proven reliable 
when applied to the diagnosis of five patients, for routine use in 
practice {i.e., as an adjuvant clinical diagnostic standard) the 
accuracy should be confirmed statistically in a larger patient 
cohort. The present study population included only 45 patients, 
and it is expected that through cooperative efforts between 
hospitals, large amounts of data can be obtained that will allow 
verification of the reliability of this method. 

Selection of EEG signal features. As many EEG signal features can 
be used as EEG parameters, the selection of suitable EEG 
characteristics as parameters is also important for the SESA 
method. A good EEG parameter represents brain injury states well 
and has values that can be easily calculated. The ApEn and a defined 
EEG frequency domain (SWC) were chosen as EEG parameters. 
Characteristics of nonlinear EEG dynamics include ApEn, sample 
and permutation entropies, and Lempel-Ziv complexity^^"^^. Future 
research will test different features to determine the best for various 
applications. In terms of EEG frequency domain features, the ratio of 
the slow wave is a good EEG feature for patients with severe disorders 
of consciousness. Whether SWC is a good EEG feature for patients 
with normal consciousness also requires further experimental 
investigation. 

The spatial resolution of SESA mapping. We used a 40 -channel 
EEG acquisition instrument to acquire EEG signals, using 16 
channels for the actual scalp electrodes. Applying high-density 
EEG signal recording instruments (with 128 or 256 channels) will 
greatly increase the spatial resolution of injured areas, facilitating 
studies that are both clinically relevant and provide basic 
functional information by revealing the topography of injured 
brain areas. However, signal interference from the various leads 



Table 4 


Probabilities of type 1 and 


1 errors 














Features 




Fp1/Fp2 




F7/F8 




T3/T4 


T5/T6 


Status 


Type 1 


Type II 


Type 1 


Type II 


Type 1 


Type II 


Type 1 Type II 


ApEn 


Quiet 


0.038 


0.048 


0.046 


0.017 


0.022 


0.037 


0.012 0.048 


Call Name 


0.032 


0.034 


0.019 


0.018 


0.042 


0.03 


0.032 0.024 


SWC 


Quiet 


0.021 


0.019 


0.039 


0.022 


0.023 


0.043 


0.045 0.030 




Call Name 


0.01 1 


0.036 


0.026 


0.041 


0.019 


0.039 


0.029 0.022 
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Table 5 


Characteristics and 


results for one patient 














Features 








Cp 


values of individual 


channel 








oTare 


FP1/FP2 


F7/F8 


U/ 1 4 


JR/JA 
1 3/ 1 O 


r J/ 




P'? /P/( 


\J \ 1 \JL 


ApEn 


Resting 


0.731 


0.788 


0.654 


0.681 


0.773 


0.815 


0.863 


0.917 


Call name 


0.702 


0.691 


0.624 


0.586 


0.786 


0.773 


0.843 


0.902 


SWC 


Resting 


1.274 


1.491 


1.657 


1.674 


1.164 


1.31 1 


1.403 


1.197 




Call name 


1.516 


1.453 


1.717 


1.786 


1.392 


1.503 


1.461 


1.265 



presents a technical limitation that must be addressed in order to 
maximize the benefits of improved spatial resolution. 

In conclusion, we innovated a SESA method to determine the 
injured areas of the brain. Future studies will focus on the develop- 
ment of specific data analysis software and additional clinical veri- 
fication, which will allow SESA to be used as an adjuvant brain 
mapping technique along with conventional imaging methods. 

Methods 

All experiments were carried out in accordance to the approved by the Wu Jing 
Hospital Institutional Review Board in Hangzhou, Zhejiang, China; and informed 
written consent was obtained from the families of all patients. 

Subjects. Subjects in this paper were placed in patient or control groups. All patients 
had unilateral brain injuries with severe disorders of consciousness, including 
minimally conscious and vegetative states^^. Patients with brain injury but who 
exhibited a normal state of consciousness were not included in the present study, but 
have been examined in other studies^" and will be the focus of future research. The 
cases in these experiments met the inclusion criteria and exclusion criteria detailed in 
Appendix I. The 45 patients were collected from inpatients received by the 
Rehabilitation Center of Zhejiang Wujing Hospital from September 201 1 to October 
2012. There were 25 cases with left brain trauma and 20 cases with right brain trauma 
(detailed information is listed in Appendix II). The control group consists of 10 
volunteers (6 males and 4 females ranging in age from 20 to 24 years of age) from 
Hangzhou Dianzi University. All of the control patients were healthy, with no mental 
disease or brain trauma. The average time between injury and study participation is at 
least 3 months; the range of time is from 3 months to 6 months. 

EEG acquisition and noise reduction. The EEC acquisition instrument with 40 
channels (Nuo Cheng, Shanghai, China) had 16 channels used to record EEG signals. 
Scalp electrodes were placed according to the international standard 10/20 system^ 
and the electrode impedance was <5 KQ, and the sampling frequency was 256 HZ, 
Al and A2 were the reference electrodes and Fpz was the grounding electrode. 

• EEG signal acquisition. EEG signals in two brain states were recorded, one under 
the resting state and the other under the call-name stimulus state. The two states were 
recorded to determine which brain state's EEG signals were better for extracting 
features for the ratio calculation (further details in the Discussion section). 

For recording in the resting state, 10 min sustainable EEG signals were acquired 
from each subject. For recording the in call-name stimulus state^^ the acquisition 
process is shown in Figure 5. There was a 12 s call- name stimulus recording, during 
which a selected patient's family member would call the patient's name three con- 
secutive times. When the third call ended, data collecting continued for 2 min before 
ending. 



• EEG noise reduction. Unstable data caused by external interference (such as EMG 
interference caused by sudden movements of patient limbs or EOG interference when 
patients blink during the acquisition process) was manually removed. The toolbox- 
IIR filter in EEGLAB was processed to remove interference caused by the 50 Hz 
power frequency. 

Extraction of EEG features. Two EEG features were chosen as parameters to 
compute ratios of symmetrical channels for injured and uninjured brain areas (and 
channels). One ratio is Approximate Entropy (ApEn)^*, and the other is a defined 
EEG feature termed the slow wave coefficient (SWC). The extraction processes of the 
two features are described below. 

• Slow wave coefficient (SWC) calculation. EEG rhythm is considered an overall 
indicator of the excitability of the central nervous system. In general, a lower 
frequency of slow waves corresponds to a greater degree of abnormality in the 
corresponding areas^^. thus, brain damage causes an increase in the slow wave and 
decrease in the fast wave band in the EEG power spectrum^'^.As such, SWC provides a 
good index for the identification of damaged brain areas. 

Since the EEG power spectrum is typically divided into six bands, i.e., 6 (1.0- 
4.0 Hz), G (4.1-8.0 Hz), oci (8.1-10.0 Hz), 0^2 (10.1-13.0 Hz), Pi (13.1-17.5 Hz), and 
(32 (17.6-30.0 Hz), SWC was defined as: 



SWC = 



spectrum of(6 + 0) 
spectrum of(a+ P) 



In this formula, a = oci + 0^2, P = Pi + P2. 
Calculation of SWC was as follows: 

(1) 2 min EEG signals of each subject in the resting state and 12 s in call-name 
stimulus state were intercepted; 

(2) The power spectrum of each band was extracted by Fourier transform; 

(3) By using the above-mentioned definition of SWC, values of each channel for 
every subject were obtained. 



• Approximate Entropy (ApEn) calculation. Approximate Entropy (ApEn) is a 
nonlinear dynamic character used to describe the complexity and regularity of signals 
to measure the irregularity of a time series^^'^*^. Its physical significance is the 
probability of generating a new model of a time series when the dimension changes. 
The larger probability indicates the sequence is more complex, leading to a greater 
ApEn. Recent research has shown that nonlinear analysis can directly monitor and 
measure the suppression level of the cerebral cortex in real-time. For patients with 
severe disorders of consciousness, the higher entropy indicates the patient is more 
conscious; thus ApEn is another index that reflects damaged brain areas. 

For the calculation of ApEn a time window was put on each channel of EEG data 
and 2 s (N = 512) chosen as the window time, then one sampling point was used as a 
window shift to calculate the ApEn of each channel. After that, the average ApEn 
from all time windows was taken as the ApEn value. Details of the computing method 
are provided elsewhere^'*. 




Figure 4 | Injured channels determined by the SESA method (left) and CT images of injured area (right) for one patient. 
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Figure 5 | Acquisition procedure under the call-name stimulus. 



Grouping of symmetrical channels. For the grouping of symmetric channels, the 
side and area of injury in each patient were identified from CT or MRI images by a 
clinician. Channels corresponding to injured areas were referred to as injured 
channels (with the offside channel designated as the symmetric uninjured channel; 
see Fig. 3). Channels for the whole brain were then divided into injured and uninjured 
areas (shown as the areas above and below the red line in Fig. 3, respectively), with 
ratios values of symmetric channels calculated as [injured channels vs. symmetric 
uninjured channels] and [uninjured channels vs. symmetric uninjured channels], 
respectively, (as shown the area under the red line in Figure 6). 

Computation of ratio values. As the ratio values are critical for the SESA method 
(according to our proposed hypothesis) the ratio values of EEC features between 
injured and uninjured brain areas as well as channels were calculated. For ApEn and 
SWC, the calculations were as follows: for injured areas, the feature (ApEn or SWC) of 
the injured channel was divided by the feature of the symmetric uninjured channel; 
for uninjured areas, the feature of the uninjured channel was divided by the feature of 
the symmetric uninjured channel. 

Statistical analysis of ratio values in population and individual channels 

(1) After ratio values were acquired, differences in ratio values for population 
{i.e., total injured vs. uninjured channels) and individual injured channels 
{i.e., injured individual Fl vs. total uninjured channels) were evaluated. For 
convenience, Cp was used instead of ratio values. The analysis of Cp values in 
population channels was performed as follows: 

(2) Computation of the mean ratio for all injured area channels; 

(3) Computation of the mean ratio for all uninjured channels; 

(4) Placement of 25 means (cases) in a column and one-way ANOVA of the ratio 
between injured and uninjured channels. 

The analysis of Cp values in individual injured channels was performed as follows: 

(1) The totals for individual injured channels were calculated as 20 FPl or FP2, 
18 F7 or F8, 17 T3 or T4, 5 T5 or T6; 

(2) The mean ratio for all individual injured channels was calculated; 

(3) The mean ratio for all individual uninjured channels was calculated; 

(4) One-way ANOVA was performed for injured and uninjured channel groups. 
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Figure 6 | Diagram of channel grouping. Two kinds of symmetrical 
channel pairs are shown, corresponding to injured (i.e., injured channel vs. 
symmetrical uninjured channel; broken dark yellow line) and uninjured 
{i.e.y uninjured channel vs. symmetrical uninjured channel; broken green 
line) areas. 



Establishing criteria for diagnosing an injured area. If the statistical analysis of ratio 
values revealed significant differences in population and individual channels between 
injured and uninjured brain areas, a determination of whether a brain injury had 
occurred was made based on a range of ratio values. 

SPASS analyses were first used to determine whether the Cp values of all channels 
follow a normal distribution. By applying statistical theory, preliminary criteria for 
judgment were derived. For an injured channel (e.g., Fpl or Fp2), the ratio of 20 
patients were taken as sample values to estimate the range of the ratio for corres- 
ponding injured brain areas (for other injured channels, the computing methods were 
the same). For an uninjured channel (e.g., F3-F4, C3-C4, P3-P4, 01-02), the ratios 
of 25 patients were regarded as sample values to estimate the range of the ratio for 
corresponding uninjured brain areas. More details about the computing methods are 
in Appendix III. 
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